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We use non-perturbative renormalization group techniques to calculate the momentum depen- 
dence of thermal fluctuations of graphene, based on a self-consistent calculation of the momentum 
dependent elastic constants of a tethered membrane. We find a sharp crossover from the perturbative 
to the anomalous regime, in excellent agreement with Monte Carlo results for the the out-of-plane 
fluctuations of graphene, and give an accurate value for the crossover scale. Our work strongly sup- 
ports the notion that graphene is well described as a tethered membrane. Ripples emerge naturally 
from our analysis. 
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I. INTRODUCTION 

Free standing graphene, the only presently known 
mono-atomic two-dimensional crystal"'^, should be an 
ideal tethered membrane, i.e. a membrane made up 
of constituents with a fixed connectivity that give rise 
to a finite shear modulus of the membrane. Tethered 
membranes are known to have highly unusual proper- 
ties, such as the absence of any finite elastic constants 
in the thermodynamic limit, a negative Poisson ratio, 
and fluctuations characterized by a large anomalous di- 
mension rj in the infra-red (IR) limilj^l. Experiments 
have not yet probed the large wavelength regime to test 
these predictions, however a negative Poisson ratio and 
anomalous fluctuations have been seen in Monte Carlo 
(MC) simulations of graphene based on an realistic ef- 
fective many-body interaction of C-atom^. One of the 
most surprising outcome of experimental investigations 
of free standing graphene were the observation of ripples 
in graphene sheets with a characteristic scale 50—100 A^. 
While it was often argued that these ripples are not 
compatible with the standard continuum elastic theory 
of tethered membraned^, we demonstrate below that 
they emerge naturally from it. Since ripples are a fi- 
nite scale phenomenon, this requires to go beyond the 
asymptotic regime which was investigated in previous 
theoretical investigations. Here, we present the most 
thorough RG treatment of tethered membranes yet, a 
non-perturbative renormalization group (NPRG) analy- 
sis of tethered membranes which is based on the expan- 
sion of the effective action in terms of elastic coupling 
functions which for the first time allows to extract the 
full momentum dependence of thermal fluctuations. Ex- 
cellent agreement with MC simulations of free standing 
graphene is found. Ripples emerge as the real space ana- 
log of the Ginzburg scale, which is the crossover scale 
which seperates the anomalous regime from the perturba- 
tive one. We further calculate the anomalous dimension 
both in the flat phase and at the crumpling transition 
which is found to be of second order. 

In contrast to fluid membranes which are always crum- 
pled (the average of the normal of the surface vanishes) 
the finite shear modulus stabilizes a flat phase with long 



range order of the normalJ^ and the normal-normal cor- 
relation function Gjq decays asymptotically for small mo- 
menta q as G'jv(g) ~ 9"^+'', see Refs. [2:8 9J. The flat 
phase is stable for ry > and in fact all calculations yield a 
large anomalous dimension, varying between rj — 2/d ~ 
2/3 from a large d expansion for Z3-dimensional mem- 
branes embedded in d-dimensional spac^, rj » .821 
from the self-consistent screening approximatiorP^l^iJ and 
rj « 0.85 where the last result was obtained both from 
NPRGflS as well as from MC simulations of graphen^. 
The analysis from Ref. [R] was able to reproduce all 
previously known results for general d,D obtained via 
perturbative RG methods within a unified framework. 
While the analysis in Ref. [H] is restricted to the lead- 
ing order of a derivative expansion of the action, here we 
will significantly extend their analysis in a way which al- 
lows to investigate thermal fluctuations at all momenta 
up to the ultraviolet (UV) cutoff (we keep here only to 
the physical most relevant situation d = D + I but our 
results are easily extended to the general case) . Our anal- 
ysis also yields the correlation functions for the in-plane 
modes, but we confine the discussion to Ghh which for 
small q is related to the normal-normal correlation func- 
tion via Gn — q^Ghh^ However, since Ghh is more read- 
ily accessible than Gat in the NPRG approach, we will 
base our analysis on Ghh- 



II. THE MODEL 

We start from a Landau-Ginzburg type ansatJ^EISEll 
for the energy functional of a tethered membrane 'H = 
+ "W^ which consists of a bending part 

n''=^J d^X {dadaRf (la) 

and a stretching part 

-H^' = j d^xf-^{daRf + \{daR-d,Rf 

+ \{daR-daRf], (lb) 
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where i? is a D + 1 dimensional vector parametrizing the 
D dimensional membrane which is embedded in a -D + 1 
dimensional space. The presence of an UV cutoff Ap is 
implicitly assumed. The inverse temperature (3 — l/ksT 
is absorbed in the definition of the effective parameters, 
i.e. K = /3k, /2 — P^, A — /3X and fo — Pfo- If one 
writes the stretching part of the membrane in terms of 
derivatives, rUa = daR{x), a — 1 ... D, the analogy to a 
Ginzburg-Landau expansion becomes apparent and one 
would expect a phase transition near f o ~ from a sym- 
metric, crumpled phase with (m^^ — which exists for 
positive fo to a symmetry broken flat phase, character- 
ized by the order parameter (m-a) = mo o = Js^a 7^ 
where J is the magnitude of the order parameter and Bq, 
a = 1 . . . D, are unit vectors which span the membrane, 

(R) = JXaBa- 

Here we shall be interested in the flat phase and there- 



fore rewrite Eq. ( lb ) by introducing the flat metric tensor 



9ab = "^0:0 • mfc^o = J^Sab- Defining gab ^ rUa ■ rUb and 

Uab = {gab - 9ab)l'^^ "^^ find, up to constanP, 



H"* = / d 



(2) 



where we used the mean field result Jaq = [— fo/(/i + 
r'A/2)]^/^ for the order parameter to cancel a term lin- 
ear in Uaa- It is convenient to separate in-planc and 
out-of-plane deformations of the membrane by intro- 
ducing Amo = daR — mo,o with Amo = {daU, dah) 
such that h corresponds to out-of-plane deformations 
and u = u°'ea. In these variables we have Uab = 
(l/2)( J9ow'' + Jdbu'' + dahdbh + daU ■ dbu). Note that if 
one ignores terms of third and fourth order in u and keeps 
derivati ves oi u only up to second order, one finds, after a 
rescalin^^'SlII! (h^u) J{h,u), the minimal model for a 
flat membrancP, U ^ IJ d^x[R{dadah)'^+flul^ + ^ula\ 
with Uab ~ dau'' + dbU°- + dahdbh and ^ — /i J"*, A — A J** 
and K = kJ^ . The minimal model does however not pos- 
sess the full symmetry of the Ginzburg-Landau model 
defined by Eqs. (la lb ) and cannot describe the crum- 



pling transition. Furthermore, neglecting the fourth or- 
der derivative terms of u prevents an accurate descrip- 
tion of the membrane fluctuations at finite momenta. We 
therefore do not use the minimal model here. 



function i?A removes IR divergences arising from modes 
with fc < A and will be specified below. The NPRG anal- 
ysis of Ref. [12] was restricted to the flow of the parame- 
ters which appear already in . While this is sufficient 
to discuss the asymptotic regime at vanishingiy small mo- 
menta, the RG transformation will in general lead to a 
momentum dependence of k, /i, and A which must be ac- 
counted for in a proper analysis of thermal fluctuations 
at finite momenta. For A < Ao we shall therefore make 
a non-local ansatz of the form Fa 



Vf with 



d^x d^x'kAix - x')dlR{x)dlR{x') (4a) 



and 



Tf = / d^x d^'x' 



MA (a; - x')Uab{x)Uab(x') 



1 r 



AA(a; - x')Uaa{x)Ubb{x') 



(4b) 



This rather natural generalization of the effective action 
allows to account for non-local correlations but at the 
same time ensures that, as long as only the coupling func- 
tions KA, MA and Aa and the parameter Ja are renormal- 
ized, the effective action retains at all A the full symmetry 
of the original model and thus all Ward identities will be 
obeyed. Apart from the approximation that we only take 
into account the irreducible correlations explicitly defined 
through Eqs. ( 4a|4b ) , which uniquely fixes the RG flow 
equations, no further approximations will be made. 

^-\r\j^ = 'x/x.^'Xj^ + 1/2 o/X/^^ + 1/2 o/v^^^ = 
'Vrv^'X/X. = S +1/2 + 1/2 



III. NON-PERTURBATIVE RG APPROACH 

The NPRG is based on the exact flow equatior(Sl for 
the cutoff dependent effective action Fa which for A = Ao 
coincides with the bare action 



^Fa 
dA 



-Tr 



a^FA 



Ra 



dRA 
dA 



(3) 



where A is the running IR cutoff and (j), (j)' are any of the 
fields or h. The trace stands for an integral over mo- 
mentum and a sum over internal indices. The regulator 



FIG. 1: (Color online) Flow diagrams for the one point vertex 
(1st Hne) and the self-energies Tj^h (2nd hne) and Eat. Forcing 
the one-point vertex to vanish for all A yields the flow of the 
order parameter JA.'^Wiggly lines correspond to Gab prop- 
agators and solid Unes to Ghh- Lines with a dash correspond 
to single-scale propagators defined via Gab = —G^bdhRh and 
Ghh = —GhhdARA- Small open circles denote the order pa- 
rameter and the solid dot a derivative with respect to A. 



To derive the NPRG equations we expand Fa in the 
fields Amo — {daU,dah). The Dyson equation for 
the Greens function Ghh of the h field, defined via 
(hqh^qf) = V6q^q'Ghh{q) whcre V is the Z?-dimensional 
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volume, is 

Gll{q) = G-\{q) + ^hh{q) (5) 

with (here and below we suppress in our notation the A 
dependence of the coupling parameters) 

'^hh[q) = (k, -KAo)g'*, (6a) 
G^X = ^Ao'z' + i?A(<z), (6b) 

where denotes the bare and momentum independent 
value of the initial coupling constant k defined at the UV 
cutoff Aq. The Greens functions of the in-plane modes, 
defined via (wfcW^Lfc/) = y^k.k'Gabik)^ can be written in 



terms of transverse and longitudinal components, 

Gab{k) = GAk) {6ab - k^k'/k^) + G\\{k)k'^k'' 1^ (7) 
with G'^^ = Gq + Sq for a =_L, ||, and 

Ei(fc) = J^fikk^ + {kk-iiAo)k\ (8a) 
S|,(fc) - J^2ilk + ~Xk)k^ + ikk-k^,)k\ (8b) 

To determine the flow of J\ and the self-energies, we 
further need the three and four point vertices. In a sym- 
metrized form they read 



-'ac ) 



rifei(^l)^2,fe3) = -ij\>^kAk\ ■ k2)k^Sab + Afc2(fci • k3)k26ac + Afe, (fc2 ' k3)k16bc + (fci • k2){flkik^Si,c + jlk^klSa 

+{ki ■ k3){ilkik2Sbc + /ife3^2'5a6) + (^2 ' k:i){ilk^k\5ac + ilkiklSab)^, (9a) 
ril(gi,g2;fc) = ~ij{fik[{qi-k)q'i + {q2-k)ql\+\k{qi-q2)k^]. (9b) 

^^abcdiki,k2,k^,kA) = M/ci2[(fcl • fc3)(fc2 ' ^4) + (^1 ' fc4)(fc2 • fe3)]^ab(5cd + Afci., [(^1 ' fe2)(fe3 ' k^) + (fci • A;4)(A;2 • k:i)]5ac5bd 
+flkiA{kl ■ k2){k3 ■ ki) + (fei • fc3)(fc2 • k4)]6adSbc 

+Xki2{ki ■ k2){k3 ■ k4)6abScd + Afei3(fci • k3){k2 ■ k4)SacSbd + Afei4(fci • k4){k2 ■ ksjSadSbc, (9c) 

ri'lab(9i'92;fei,fe2) = '^a6{A9i2 [(gi • fei)(92 ■ fe2) + (gi ■ fc2)(g2 • fei)] + ■ g2)(fci • fe2)} , (9d) 
92, 93, 94) = (Agi2 + + Agi3)(gi • '73)(g2 • qa) + {h,^ + Agia + \i4)(9i • 94)('72 • 93) 

+(Mgi3 + /^9i4 + Agi2)(9i • 92)(g3 • 94), (9e) 



where fcy = |fci + The subscript /i and momenta q.; 
refer to h fields while subscripts a . . .d and fci refer to tt 
fields. The flow equations for the order parameter J and 
the self-energies T,h,h, and Em are rather long and 
shown diagrammatically in Fig. [ij They yield coupled 
integro-differcntial equations for the flow of the coupling 
functions which must be solved self-consistently. Note 
that the equations are closed, since all three- and four- 
point vertices in Eqs. ( 9a]|9e | are entirely determined by 
coupling functions which can be extracted from the self- 
energies. This is a result of the form of Fa in Eqs. ( 4a|4b ) 
which relates the third and fourth order vertices to lower 
order ones which in a usual field expansion would have 
to be imposed through Ward identities. 

We integrate the flow equations from the UV cutoff 
A = Aq to A = numerically, using for numerical stabil- 
ity an analytic regulator, RA{q) = Kq=o q^ I (exp[(g'/A)^] — 
1) (this one, as well as a non-analytic regulator, were 
also used in^-'). If we ignore the momentum dependence 
of Kg, liq, and \q our flow equations reduce to those re- 
ported in Ref. I_l2i|. In the flat phase we find for Z? = 2 
r\ « 0.85 which agrees with the derivative expansion 
resultP. For completeness, we note that our NPRG ap- 



proach yields for D = 2 a second order crumpling tran- 
sition (to within numerical accuracy) with an anomalous 
dimension 77 k, 0.64(5), slightly larger than the result 
T] K, 0.627 obtained with a sharp cutoff and a derivative 
expansion in^^, where a weak dependence of the flow on 
the form of the regulator prevented a firm conclusion on 
the order of the transition. Our result, obtained with 
a smooth regulator, resolves this ambiguity in favor of 
a transition of second order. However, we cannot rule 
out that terms of higher order in the stress tensor would 
change the nature of the transition. 



IV. COMPARISON WITH MC DATA AND THE 
ROLE OF THE GINZBURG SCALE 

The MC data for Ghh were obtained from a system 
of 37888 C-atoms with an accurate b ond ord er potential 
LCBOPII and T = 300 K, see Refs.l2lSIMi3 for details. 
To reduce statistical noise, the out-of-plane distortions 
hi were obtained by evaluating for each i the average 
hi = (3/ii -I- ^i)/6 where the sum runs over the 

three neighbors of atom z.l^ As expected, for small q one 
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finds the relation q^Ghh — Gn between the correlation 
function of the normals and the height fluctuation which 

for graphene is extremely accurate even up to q w 1 A , 
see the inset in Fig. [2] Since the very small q data for 
Ghh is more noisy than that of Gat we used for the last 
three data points (q < 0.07 A ) the data of Gn to calcu- 
late Ghh, the result is shown in Fig. |2] The strong peak 
near = 47r/3a ~ 2.94 A , where a is the equilibrium 
lattice parameter, corresponds to the first Bragg peak. It 
defines the upper limit beyond which continuum theory 
is inapplicable and it serves as a natural UV cutoff for 
the NPRG calculation, Ag = gs. For smaller q the data 
shows the scaling Ghh oc q~'^ of the perturbative regime 
and for very small q the anomalous scaling Ghh oc (7""*+'' 
where 77 ~ 0.85 agrees with the NPRG result. 

An important scale is the Ginzburg scale qQ for the 
crossover from the perturbative to the anomalous regime. 
This scale cannot be captured within a finite order deriva- 
tive expansion. Perturbation theory^ yields for D = 2 
the rough estimate 



pt 



1g 



[3ifo/(27r)]i/V(4^) 



(10) 



with Kq = AJ + \)/{2ji + A). Below we extract 
qQ from the numerical flow which gives a more accurate 
estimate. 




FIG. 2: (Color online) Results for the out-of-plane fluctua- 
tions q^Ghh{q) from NPRG (solid, black), MC (dashed, red) 
and from a simple phenomenological approximation discussed 
in the main text (dashed-dotted, blue). The inset shows MC 
data for the out-of-plane correlations q^Ghhil) (dashed, red) 
and the fluctuations of normals Gn^q) (solid, green). 



The RG equations require the initial form of kg, fig, 
and the mean-field order parameter Jaq which are de- 
fined at the UV cutoff Aq. While we cannot rule out a 
q dependence of the initial coupling functions, for sim- 
plicity and in accordance with Eqs. (la lb I we choose 



g-independent constant values for the initial form of the 
coupling functions. To fix the initial value of k we us e a 
value previously reported in the literature, k = 1.1 eV.!^ 
The elastic properties were studied in detail in Ref. [16] 
and values for the bulk modulus B = J^(p + X) and 



the shear modulus, in our notation J^/i, were extracted 
for moderate system sizes. For scales smaller than the 
Ginzburg scale, all elastic constants are strongly cutoff 
dependent and in particular in the IR limit Aq=o and 
f^q=o vanish as \*~D-2ri .^^j^gj-gag k^^q diverges as A"**. 
Values of the clastic constants of free standing graphene 
must therefore be understood as valid only for a given 
system size (or IR cutoff). The system size used in 
Ref. [TB] are for T = 300 K at the border of the pertur- 
bative regime and the reported elastic constants are not 
yet strongly renormalized and close to those at smaller 
scales. We therefore use these results to fix the initial val- 
ues of J^fi « 9.95 eVA"^ and .PX « 2.41eVA"^ The 
initial value of the order parameter J^^ « 2.5 is chosen 
to give the best overall agreement with the MC data. 
These values place graphene well inside the fiat phase, 
in accordance with numerical simulations which show 
no sign of crumpling even at high temperature^^. The 
NPRG result for Ghh with these values and T = 300 K 
is shown in Fig. [2j The agreement with the MC data 
is very good, especially the sharpness of the crossover 
from the perturbative to the anomalous regime is well 
reproduced. A simple phenomenological ansatz for the 
self-energy, Y^hh = fiAo^Q^iQc/Q)^ with A fixed by the 
asymptotic behavior, does not lead to a satisfactory de- 
scription of the data, as was already noted in Ref. 
see Fig. [2] The Ginzburg scale is by standard defini- 
tion the scale where selfenergy correction to kq^o equal 

the bare parameter which allows to read qc ~ 0.08 A 
directly off the flow, which is slightly smaller than the 

perturbative estimate q^ ~ 0.12 A from Eq. (10). The 
Ginzburg length Lq — 27: /qo « 80 A is of the same order 
as experimentally observed rippleP' which offers a natu- 
ral explanation of their appearance as just the real-space 
manifestation of the Ginzburg scale. Furthermore, ex- 
perimentally the fluctuations were found to be broadly 
distributed around a characteristic scale, which is again 
in accordance with the behavior of Gh around qc which 
is not characterized by a sharp feature at qo but by a 
gradual crossover from the perturbative to the anomalous 
regime. The qualitatively correct perturbative estimate 
of the Ginzburg scale Eq. ( 10 1 furthermore yields a sim- 



ple dependence of the Ginzburg scale on the temperature, 
qc ~ ^/KqUbT / K. If our interpretation of ripples in free 
standing graphene is correct, the average real space scale 
of ripples should thus increase as k/^/KqIcbT on lower- 
ing the temperature, where Kq is the two dimensional 
Young's modulus. 



V. CONCLUSIONS 

In summary, we have presented a NPRG analyis for 
tethered membranes which avoids a derivative expansion 
and is the first to include the full momentum dependence 
of the elastic coupling parameters. Our solution of the 
NPRG fiow equations is completely self-consistent and 
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obeys all symmetry constraints. In our approach the 
crumpling transtion is found to be of second order and we 
give an improved estimate for the anomalous dimension 
at the transition. For the flat phase of the membrane, 
we find excellent agreement with MC results for the mo- 
mentum dependence of the out-of-plane fluctuations of 
free standing graphene. Also the crossover region, which 
shows a relatively sharp crossover from the perturbative 
regime to the anomalous scaling regime which is char- 
acterized by a large anomalous dimension, is well repro- 
duced. This strongly supports the notion that free stand- 
ing graphene behaves just as a tethered membrane, albeit 
a very stiff one. The most important scale in the analysis 



of the momentum dependence of the membrane fluctu- 
ations is the Ginzburg scale which we find to be of the 
same order as the experimentally determined characteris- 
tic size of ripples. The observation of ripples at this scale 
should thus be looked at as a confirmation of the con- 
tinuum elastic theory of tethered membranes, a notion 
which could also be tested experimentally by measuring 
the characteristic ripple scale as a funtion of temperature. 
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